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Abstract 

The integrand of any multi-loop integral is characterised after Feynman parametrisation by 
two polynomials. In this review we summarise the properties of these polynomials. Topics 
covered in this article include among others: Spanning trees and spanning forests, the all- 
minors matrix-tree theorem, recursion relations due to contraction and deletion of edges, 
Dodgson's identity and matroids. 



1 Introduction 



In this review we discuss Feynman graph polynomials. Let us first motivate the interest in these 
polynomials. The Feynman graph polynomials are of interest from a phenomenological point 
of view as well as from a more mathematical perspective. We start with the phenomenological 
aspects: For the practitioner of perturbative loop calculations the integrand of any multi-loop 
integral is characterised after Feynman parametrisation by two polynomials. These two polyno- 
mials, which are called the first and the second Symanzik polynomial, can be read off directly 
from the Feynman graph and are the subject of this review. They have many special properties 
and these properties can be used to derive algorithms for the computation of loop integrals. 

In recent years the graph polynomials have received additional attention from a more formal 
point of view [fH-l29l. Feynman integrals are now considered as non-trivial examples of mixed 
Hodge structures and motives. The zero sets of the graph polynomials play a crucial role in this 
setting. 

Graph polynomials have a long history dating back to Kirchhoff ll30ll . There are well- 
established books on this subject II3TH361I . However the field has evolved, new insights have 
been added and old results have been re-discovered. As a consequence the available information 
is scattered over the literature. What is missing is a concise summary of the properties of these 
polynomials. With this review we hope to fill this gap. We have tried to make this article acces- 
sible both to the phenomenological oriented physicist interested in loop calculations as well as 
to the mathematician interested in the properties of Feynman integrals. In addition we include in 
a few places results which are new and cannot be found in the literature. 

This review is organised as follows: In the next section we recall some basic facts about 
multi-loop Feynman integrals. We introduce the two polynomials associated to a Feynman graph 
and give a first method for their computation. Sect.[3]is devoted to the singularities of a Feynman 
integral. There are two aspects to it. If we fix the external momenta and take them to lie in 
the Euclidean region the singularities of the Feynman integral after integration are of ultraviolet 
or infrared origin. They arise - apart from a possible overall ultraviolet divergence - from the 
regions in Feynman parameter space where one of the two graph polynomials vanishes. For the 
second aspect we give up the restriction on the Euclidean region and view the Feynman integral 
as a function of the external momenta. As we vary the external momenta, additional threshold 
singularities may arise. Necessary conditions for them are given by Landau's equations. In 
sect. ID we start to introduce concepts of graph theory and define for a graph its spanning trees 
and spanning forests. These concepts lead to a second method for the computation of the graph 
polynomials, such that the two graph polynomials can be directly read off from the topology of 
the graph. Sect.[5]introduces the Laplacian of a graph and states the matrix-tree theorem. This in 
turn provides a third method for the computation of the two graph polynomials. This method is 
well suited for computer algebra, as it involves just the computation of a determinant of a matrix. 
The matrix is easily constructed from the data defining the graph. In sect. [6] the two operations 
of deleting and contracting an edge are studied in detail. This leads to a fourth and recursive 
method for the computation of the two graph polynomials. In addition we discuss in this section 
the multivariate Tutte polynomial and Dodgson's identity. Sect. [7] is devoted to the dual of a 
(planar) graph. The Kirchhoff polynomial and the first Symanzik polynomial exchange their role 
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when going from a graph to its dual. Sect.[8]is of a more formal character. We introduce matroids, 
which provide a natural generalisation of graphs. Within matroid theory, some things are simpler: 
For example there is always the dual of a matroid, whereas for graphs we were restricted to planar 
graphs. Matroid theory provides in addition an answer to the question under which conditions 
two topologically different graphs have the same Kirchhoff polynomial. Finally, sect. |9] contains 
our conclusions. 

2 Feynman integrals 

In this section we recall some basic facts about multi-loop Feynman integrals. We introduce the 
two polynomials associated to a Feynman graph and give a first method for their computation. 
We will work in a space-time of D dimensions. To set the scene let us consider a scalar Feynman 
graph G with m external lines and n internal lines. Fig. Q] shows an example. In this example 
there are four external lines and seven internal lines. The momenta flowing in or out through the 
external lines are labelled p\, p 2 , pi and p4 and can be taken as fixed D-dimensional vectors. 
They are constrained by momentum conservation: If all momenta are taken to flow outwards, 
momentum conservation requires that 

PI+P2+P3+P4 = 0. (1) 

At each vertex of a graph we have again momentum conservation: The sum of all momenta 
flowing into the vertex equals the sum of all momenta flowing out of the vertex. 

A graph, where the external momenta determine uniquely all internal momenta is called a 
tree graph. It can be shown that such a graph does not contain any closed circuit. In contrast, 
graphs which do contain one or more closed circuits are called loop graphs. If we have to specify 
besides the external momenta in addition I internal momenta in order to determine uniquely all 
internal momenta we say that the graph contains / loops. In this sense, a tree graph is a graph with 
zero loops and the graph in fig.[T]contains two loops. In more mathematical terms the number / is 
known as the cyclomatic number or the first Betti number of the graph. Let us agree that we label 
the / additional internal momenta by k\ to k\. In the example of fig. Q] there are two independent 
loop momenta. We can choose them to be k\ = qi and k 2 = 16- Then all other internal momenta 
are expressed in terms of k\, k 2 and the external momenta p\, p^. 

q\=k\-p\, q 2 = h -pi -pi, q4 = h+k 2 , ^ 

<?5 = k 2 ~P3-pA, qi = k 2 -pA- 

In general, each momentum flowing through an internal line is given as a linear combination of 
the external momenta p and the loop momenta k with coefficients —1, or 1: 

l m 

qi = J^hjkj+J^GijPj, E {-1,0,1}. (3) 

7=1 7=1 

We associate to a Feynman graph the Feynman integral 

J r=\ 7=1 t ^7 +m 7^ 7=1 
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Figure 1 : The "double box"-graph: A two-loop Feynman diagram with four external lines and 
seven internal lines. The momenta flowing out along the external lines are labelled pi, p4, the 
momenta flowing through the internal lines are labelled q\, qj. 

This Feynman integral depends on the Feynman graph G and the external momenta p. The 
graph G together with the independent loop momenta k and the external momenta p fixes all 
internal momenta according to eq. ©. In addition Iq depends on the masses nij of the particles 
corresponding to the internal line j, as well as integer numbers Vj, specifying the power to which 
each propagator is raised. The parameter p is an arbitrary mass scale. We have multiplied the 
integral with a factor (p 2 ) v ~ lD / 2 , this ensures that Iq is dimensionless. 

How to perform the integration over the loop momenta? The first step is to convert the 
products of propagators into a sum. This can be done with the Feynman parameter technique. In 
its full generality it is also applicable to cases, where each factor in the denominator is raised to 
some power Vj. The formula reads: 

n4- = — v=2> y . (5) 




We use this formula with Pj = —q 2 + m 2 . Applied to eq. © we have 

j^XiPi = j^xii-qf + m 2 ). (6) 
;=l i=l 

Now one can use translational invariance of the D-dimensional loop integrals and shift each loop 
momentum k r to complete the square, such that the integrand depends only on k 2 . Then all D- 
dimensional loop integrals can be performed. As the integrals over the Feynman parameters still 
remain, this allows us to treat the D-dimensional loop integrals for Feynman parameter integrals. 
One arrives at the following Feynman parameter integral ll37l : 

nr(v ; ) Xj io '=! v=i J 7 
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The functions 11 and J depend on the Feynman parameters Xj. If one expresses 

£x 7 (-^ 2 + m J 2 ) = -££* r M„* J + £2* r .0 r + /, (8) 

j=\ r=ls=l r=l 

where M is a I x / matrix with scalar entries and 2 is a /-vector with four-vectors as entries, one 
obtains 

11 = det(M), 7 = det(Af) (7 + QM~ l Q) /> 2 . (9) 

The functions 11 and f are called graph polynomials and are the subject of this review. They are 
polynomials in the Feynman parameters and - as we will show later - can be derived from the 
topology of the underlying graph. The polynomials 11 and f have the following properties: 

• They are homogeneous in the Feynman parameters, 11 is of degree Z, f is of degree l + l. 

• 11 is linear in each Feynman parameter. If all internal masses are zero, then also f is linear 
in each Feynman parameter. 

• In expanded form each monomial of 11 has coefficient +1. 

We call 11 the first Symanzik polynomial and J the second Symanzik polynomial. Eqs. © and 
© allow us to calculate these polynomials for a given graph. We will learn several alternative 
ways to determine these polynomials later, but for the moment it is instructive to go through this 
exercise for the graph of fig. CD We will consider the case 

pj = 0, pl = 0, pl = 0, p% = 0, 

m l — m 2 — m 3 = m 4 — m 5 — m 6 = m l = 0- (10) 



We define 



We have 



(Pl+P2) 2 = (P3+P4) 2 , t = (p 2 +P3) 2 = (PI+P4) 2 ■ (11) 
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^Xj{— qj) = - (x\ +X2+X2, +X4)k\ —2x4k\ -k 2 - (^4+^5 +x 6 +x-])kl (12) 

+2 [xipi +x 2 (pi + pi)} -h+2 [x 5 (p3 + pa) +x-jp4} ■ k 2 - (x 2 +x 5 )s. 
In comparing with eq. ® we find 



M = 
Q = 



^1+^2+^3+^4 X4 

X4 X4 + X5 + X(, + Xj 

xipi +x 2 (pi +p 2 ) 

X5 (P3+P4) +XJP4 

J = ( X2+Xs )(- S ). (13) 
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Plugging this into eq. © we obtain the graph polynomials as 




(14) 



We see in this example that 11 is of degree 2 and J is of degree 3. Each polynomial is linear in 
each Feynman parameter. Furthermore, when we write 11 in expanded form 



each term has coefficient + 1 . 

3 Singularities 

In this section we briefly discuss singularities of Feynman integrals. There are two aspects to 
it. First we fix the external momenta and take them to lie in the Euclidean region. We may 
encounter singularities in the Feynman integral. These singularities are of ultraviolet or infrared 
origin and require regularisation. We briefly discuss how they are related to the vanishing of the 
two graph polynomials 11 and f . For the second aspect we consider the Feynman integrals as a 
function of the external momenta, which are now allowed to lie in the physical region. Landau's 
equations give a necessary condition for a singularity to occur in the Feynman integral as we 
vary the external momenta. 

It often occurs that the Feynman integral as given in eq. © or in eq. © is an ill-defined and 
divergent expression when considered in D = 4 dimensions. These divergences are related to 
ultraviolet or infrared singularities. Dimensional regularisation is usually employed to regulate 
these divergences. Within dimensional regularisation one considers the Feynman integral in 
D = 4 — 2e dimensions. Going away from the integer value D = 4 regularises the integral. In D = 

4 — 2e dimensions the Feynman integral has a Laurent expansion in the parameter e. The poles 
of the Laurent series correspond to the original divergences of the integral in four dimensions. 

From the Feynman parameter integral in eq. © we see that there are three possibilities how 
poles in e can arise: First of all the Gamma-function T(v — ID/ 2) of the prefactor can give rise 
to a (single) pole if the argument of this function is close to zero or to a negative integer value. 
This divergence is called the overall ultraviolet divergence. 

Secondly, we consider the polynomial 11. Depending on the exponent v — (/ + l)D/2 of 
11 the vanishing of the polynomial 11 in some part of the integration region can lead to poles 
in e after integration. As mentioned in the previous section, each term of the expanded form 
of the polynomial 11 has coefficient + 1 , therefore 11 can only vanish if some of the Feynman 
parameters are equal to zero. In other words, 11 is non-zero (and positive) inside the integration 



11 = X\X5 + X[X(, +XIX-J +X2X5 +X2X(, +X2X7 + X3X5 +X3X6 

+X1X4 + X2X4 + X3X4 + X4X5 + JC4X6 + X4X1 , 



(15) 
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region, but may vanish on the boundary of the integration region. Poles in e resulting from the 
vanishing of 11 are related to ultraviolet sub-divergences. 

Thirdly, we consider the polynomial J . In an analytic calculation one often considers the 
Feynman integral in the Euclidean region. The Euclidean region is defined as the region, where 
all invariants (p^ + pi 2 + ... + Pi k ) 2 are negative or zero, and all internal masses are positive 
or zero. The result in the physical region is then obtained by analytic continuation. It can be 
shown that in the Euclidean region the polynomial J is also non-zero (and positive) inside the 
integration region. Therefore under the assumption that the external kinematics is within the 
Euclidean region the polynomial J can only vanish on the boundary of the integration region, 
similar to what has been observed for the the polynomial U . Depending on the exponent v — lD/2 
of J the vanishing of the polynomial J on the boundary of the integration region may lead to 
poles in £ after integration. These poles are related to infrared divergences. 

The Feynman integral Iq as given in eq. © depends through the polynomial J on the external 
momenta pj. We can also discuss Iq as a function of the p/s without restricting the external 
kinematics to the Euclidean region. Doing so, the region where the polynomial J vanishes is 
no longer restricted to the boundary of the Feynman parameter integration region and we may 
encounter zeros of the polynomial J inside the integration region. The vanishing of J may in 
turn result in divergences after integration. These singularities are called Landau singularities. 
Necessary conditions for the occurrence of a Landau singularity are given as follows: A Landau 
singularity may occur if J = and if there exists a subset S of {1, ...,n} such that 

Xi = for i e S 
and -^-J = for j e {l,...,n}\S. (16) 

dXj 

The case corresponding to S = is called the leading Landau singularity, and cases corresponding 
to S 7^ are called non-leading singularities. It is sufficient to focus on the leading Landau 
singularity, since a non-leading singularity is the leading Landau singularity of a sub-graph of G 
obtained by contracting the propagators corresponding to the Feynman parameters x ; with i E S. 

Let us now consider the leading Landau singularity of a graph G with m external lines. We 
view the Feynman integral as a function of the external momenta and a solution of the Landau 
equations is given by a set of external momenta 

{p h p2,...,Pm} (IV) 

satisfying momentum conservation and eq. (PT6l) . If the momenta in eq. (flTT) define a one- 
dimensional sub-space, we call the Landau singularity a normal threshold. In this case all exter- 
nal momenta are collinear. If on the contrary the momenta in eq. (flTT ) define a higher-dimensional 
sub-space, we speak of an anomalous threshold. 

We give a simple example for a Feynman integral with an ultraviolet divergence and a normal 
threshold. The graph in fig. |2] shows a one-loop two-point function. The internal propagators 
correspond to particles with mass m. This graph corresponds in D — 4 — 2e dimensions to the 
Feynman integral 

(2\ E f^l \ 

~ [P } ' iiS (-k 2 +m 2 )(-(k-p) 2 +m 2 Y { } 
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Figure 2: The one-loop two-point function with equal masses. This graph shows a normal thresh- 
old for p 2 = Am 2 . 



Introducing Feynman parameters one obtains the form of eq. ©: 



Iq = r(e) / dx 



x(l —x) 



+ ■ 
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This integral is easily evaluated with standard techniques ll38l : 
Ig 
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(19) 



(20) 



Here, denotes Euler's constant. The l/e-term corresponds to an ultraviolet divergence. As a 
function of p 2 the integral has a normal threshold at p 2 = Am 2 . The normal threshold manifests 
itself as a branch point in the complex p 2 -plane. 



4 Spanning trees and spanning forests 

In this section we start to introduce concepts of graph theory. We define spanning trees and 
spanning forests. These concepts lead to a second method for the computation of the graph 
polynomials. We consider a connected graph G with m external lines and n internal lines. Let r 
be the number of vertices of the graph G. We denote the set of internal edges of the graph G by 

{ei,e 2 ,...,e n } (21) 

and the set of vertices by 

{vi,v 2 ,...,v,-}. (22) 

As before we denote by / the first Betti number of the graph (or in physics jargon: the number of 
loops). We have the relation 

/ = n-r+l. (23) 

If we would allow for disconnected graphs, the corresponding formula for the first Betti number 
would be n — r + k, where k is the number of connected components. A spanning tree for the 
graph G is a sub-graph T of G satisfying the following requirements: 
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Figure 3: The left picture shows a spanning tree for the graph of fig.ED the right picture shows a 
spanning 2-forest for the same graph. The spanning tree is obtained by deleting edges 4 and 7, 
the spanning 2-forest is obtained by deleting edges 1, 4 and 7. 



• T contains all the vertices of G, 

• the first Betti number of T is zero, 

• T is connected. 

If T is a spanning tree for G, then it can be obtained from G by deleting I edges. In general 
a given graph G has several spanning trees. We will later obtain a formula which counts the 
number of spanning trees for a given graph G. A spanning forest for the graph G is a sub-graph 
F of G satisfying just the first two requirements: 

• F contains all the vertices of G, 

• the first Betti number of F is zero. 

It is not required that a spanning forest is connected. If F has k connected components, we say 
that F is a k- forest. A spanning tree is a spanning 1 -forest. If F is a spanning fc-forest for G, then 
it can be obtained from G by deleting l + k—1 edges. Fig. [3] shows an example for a spanning 
tree and a spanning 2-forest for the graph of fig. ED 

We denote by T the set of spanning forests of G and by T k the set of spanning fc-forests of G. 
Obviously, we can write T as the disjoint union 

r 

<T = \jT k . (24) 

k=\ 

1\ is the set of spanning trees. For an element of we write 

(T u T 2 ,...,T k ) G <r k . (25) 

The T[ are the connected components of the fc-forest. They are necessarily trees. We denote by 
Pjj the set of external momenta attached to 7]. For the example of the 2-forest in the right picture 
of fig. [3] we have 

Pr 1 ={p2,P3}, Pt 2 = {puP4}- (26) 
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The spanning trees and the spanning 2-forests of a graph G are closely related to the graph 
polynomials U and f of the graph. We have 

u = e n*' 



The sum is over all spanning trees for U , and over all spanning 2-forests in the first term of the 
formula for J . Eq. (|271) provides a second method for the computation of the graph polynomials 
U and 7 . Let us first look at the formula for 11 . For each spanning tree T we take the edges 
ei, which have been removed from the graph G to obtain T . The product of the corresponding 
Feynman parameters gives a monomial. The first formula says, that U is the sum of all the 
monomials obtained from all spanning trees. The formula for J has two parts: One part is related 
to the external momenta and the other part involves the masses. The latter is rather simple and 
we write 

7 = 7o + U^xi^. (28) 

We focus on the polynomial Jq. Here the 2-forests are relevant. For each 2-forest (Ti,^) we 
consider again the edges e,, which have been removed from the graph G to obtain (T\,T2). The 
product of the corresponding Feynman parameters %i defines again a monomial, which in addition 
is multiplied by a quantity which depends on the external momenta. We define the square of the 
sum of momenta through the cut lines of {T\,T2) by 

'KM = ( I 9j) ■ (29) 

\«^(Zi,r 2 ) / 

Here we assumed for simplicity that the orientation of the momenta of the cut internal lines 
are chosen such that all cut momenta flow from T\ to T2 (or alternatively that all cut momenta 
flow from T2 to T\, but not mixed). From momentum conservation it follows that the sum of the 
momenta flowing through the cut lines out of T\ is equal to the negative of the sum of the external 
momenta of T\ . With the same reasoning the sum of the momenta flowing through the cut lines 
into T2 is equal to the sum of the external momenta of T2. Therefore we can equally write 

HtuT 2 ) = -( £ Pi) 'I £ Pj) (30) 



mr 



and jFo is given by 



^0 = i n * m. (3D 

(T h T 2 )eT 2 W(7i,r 2 ) / \ f 
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Figure 4: A two-loop two-point graph. 






Figure 5: The set of spanning trees for the two-loop two-point graph of fig.|U 

Since we have to remove / edges from G to obtain a spanning tree and (/ + 1) edges to obtain 
a spanning 2-forest, it follows that 11 and J are homogeneous in the Feynman parameters of 
degree I and (/ + 1), respectively. From the fact, that an internal edge can be removed at most 
once, it follows that 11 and Jo are linear in each Feynman parameter. Finally it is obvious from 
eq. (1271) that each monomial in the expanded form of 11 has coefficient +1. 

Let us look at an example. Fig. |4] shows the graph of a two-loop two-point integral. We take 
again all internal masses to be zero. The set of all spanning trees for this graph is shown in fig.[5J 
There are eight spanning trees. Fig. [6] shows the set of spanning 2-forests for this graph. There 
are ten spanning 2-forests. The last example in each row of fig.[6]does not contribute to the graph 
polynomial J , since the momentum sum flowing through all cut lines is zero. Therefore we have 
in this case S(t u t 2 ) = 0. In all other cases we have S(T\,t 2 ) = P 1 ■ We arrive therefore at the graph 
polynomials 



5 The matrix-tree theorem 

In this section we introduce the Laplacian of a graph. The Laplacian is a matrix constructed 
from the topology of the graph. The determinant of a minor of this matrix where the z-th row 



11 = (XI+X4)(X2+X3) + (XI+X2+X3+X4)X5, 

J = [(XI+X2)(X3 +^4)^5+^1^4(^2 +^3) +*2*3 (xi+^4)] 




(32) 
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Figure 6: The set of spanning 2-forests for the two-loop two-point graph of fig.Hl 



and column have been deleted gives us the Kirchhoff polynomial of the graph, which in turn 
upon a simple substitution leads to the first Symanzik polynomial. We then show how this 
construction generalises for the second Symanzik polynomial. This provides a third method for 
the computation of the two graph polynomials. This method is very well suited for computer 
algebra systems, as it involves just the computation of a determinant of a matrix. The matrix is 
easily constructed from the data defining the graph. 

We begin with the Kirchhoff polynomial of a graph. This polynomial is defined by 



E IT 

T(z1\ ej€T 



(33) 



The definition is very similar to the expression for the first Symanzik polynomial in eq. (1271) . 
Again we have a sum over all spanning trees, but this time we take for each spanning tree the 
monomial of the Feynman parameters corresponding to the edges which have not been removed. 
The Kirchhoff polynomial is therefore homogeneous of degree (n — /) in the Feynman param- 
eters. There is a simple relation between the Kirchhoff polynomial and the first Symanzik 
polynomial 11 : 



tl(xi,...,x n ) =x\...x n 9C 



1 

x\ 



1 

Xfi 



Qt(xi,...,x n ) =X\...X n U 



1 

XI 



1 

Xfi 



(34) 



These equations are immediately evident from the fact that 11 and %_ are homogeneous polyno- 
mials which are linear in each variable together with the fact that a monomial corresponding to a 
specific spanning tree in one polynomial contains exactly those Feynman parameters which are 
not in the corresponding monomial in the other polynomial. 

We now define the Laplacian of a graph G with n edges and r vertices as a r x r-matrix L, 
whose entries are given by H39U40H 



if i = j and edge is attached to v ; and is not a self-loop, 
H x k if i ^ j and edge connects V; and Vy. 



(35) 



The graph may contain multiple edges and self-loops. We speak of a multiple edge, if two 
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Figure 7: The left picture shows a graph with a double edge, the right picture shows a graph with 
a self-loop. 



vertices are connected by more than one edge. We speak of a self-loop if an edge starts and 
ends at the same vertex. In the physics literature a self-loop is known as a tadpole. Fig. [7] shows 
a simple example for a double edge and a self-loop. If the vertices v, and vj are connected by 
two edges and ek 2 , then the Laplacian depends only on the sum + Xk 2 . If an edge e\ is a 
self-loop attached to a vertex v,-, then it does not contribute to the Laplacian. 

Let us consider an example: The Laplacian of the two-loop two-point graph of fig.|4]is given 

by 

/ x\ +X4 —x\ 

-X\ X1+X2+X5 -X2 
— X2 X2 + Xj, 

-X4 -X 5 -X 3 



\ 



\ 



(36) 



— X4 

-X5 
x 3 +x 4 +x 5 ) 

In the sequel we will need minors of the matrix L and it is convenient to introduce the following 
notation: For arxr matrix A we denote by A[i\ , . . . , ii\ j\, . . ., jk] the (r — k) x (r — k) matrix, which 
is obtained from A by deleting the rows i\, ^ and the columns j\, jk- For A[i\, ik'Ji, h] 
we will simply write A [i\ , . . . , ik] . 

Let vt be an arbitrary vertex of G. The matrix-tree theorem states [|39l 



9C = det L[i] 



(37) 



i.e. the Kirchhoff polynomial is given by the determinant of the minor of the Laplacian, where 
the i-th row and column have been removed. One can choose for i any number between 1 and r. 

Choosing for example i = 4 in eq. (|36l ) one finds for the Kirchhoff polynomial of the two-loop 
two-point graph of fig. |4] 



(38) 



— x\ 

— X\ X1+X2+X5 —X2 
— X2 X2 + X3 

= ^1X2(^3 +X4) + (X[ +^2)^3X4 + [x\X2 + X1X3 +X2X4 +X?,X4)X5. 

Using eq. (|34|) one recovers the first Symanzik polynomial of this graph as given in eq. (|32|) . 

The matrix-tree theorem allows to determine the number of spanning trees of a given graph 
G. Setting X{ = ... = x n = 1, each monomial in and 11 reduces to 1. There is exactly one 
monomial for each spanning tree, therefore one obtains 



ac(l,...,l) = «(l,...,l). 



(39) 
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The matrix-tree theorem as in eq. (1371) relates the determinant of the minor of the Laplacian, 
where the z'-th row and the z'-th column have been deleted to a sum over the spanning trees of the 
graph. There are two generalisations we can think of: 

1 . We delete more than one row and column. 

2. We delete different rows and columns, i.e. we delete row i and column j with i ^ j. 

The all-minors matrix-tree theorem relates the determinant of the corresponding minor to a spe- 
cific sum over spanning forests ll4Tl443ll . To state this theorem we need some notation: We 
consider a graph with r vertices. Let I = (z'i,...,4) with 1 < i\ < ... < < r denote the rows, 
which we delete from the Laplacian, and let J = (ji, A) with 1 < j\ < ... < ji < r denote the 
columns to be deleted from the Laplacian. We set |/| = i\ + ... + 4 and |7| = j\ + ... + j^. We 
denote by T, ' the spanning /c-forests, such that each tree of an element of T^ ,J contains exactly 
one vertex v, a and exactly one vertex v Jp . The set l" k ' is a sub-set of all spanning /c-forests. We 

now consider an element F of ll ,J . Since the element F is a /c-forest, it consists therefore of k 
trees and we can write it as 

F = (T h ...,T k ) E ll' J . (40) 

We can label the trees such that v/, 6 T\, vt k G 7^. By assumption, each tree T a contains also 
exactly one vertex from the set {v/j, ...,Vj k }, although not necessarily in the order Vj a E T a . In 
general it will be in a different order, which we can specify by a permutation lip e S^: 

v ja eT Ma) . (41) 
The all-minors matrix-tree theorem reads then 

dot L[IJ] = £ sign(jc F ) J] *j- ( 42 ) 

-i,J e,eF 



In the special case I = J this reduces to 



detL[7] = £ n*i" ( 43 ) 



If we specialise further to I = J = (i), the sum equals the sum over all spanning trees (since 
each spanning 1-forest of q 1 ^'^ necessarily contains the vertex v ; ). We recover the classical 
matrix-tree theorem: 

det L[i] = £ n*i- (44) 

Let us illustrate the all-minors matrix-tree theorem with an example. We consider again the two- 
loop two-point graph with the labelling of the vertices as shown in fig. [8] Taking as an example 

/= (2,4) and 7= (3,4) (45) 
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Figure 8: The left picture shows the labelling of the vertices for the two-loop two-point function. 
The middle and the right picture show the two 2-forests contributing to T 2 I,J w i m I = (2,4) and 
7= (3,4). 



we find for the determinant of L[I;J] : 
det L[2,4;3,4] 



X\ + X4 —X\ 
-X2 



-X\X2~X2X\. (46) 



On the other hand there are exactly two 2-forests, such that in each 2-forest the vertices V2 and 
V3 are contained in one tree, while the vertex V4 is contained in the other tree. These two 2- 
forests are shown in fig.[8l The monomials corresponding to these two 2-trees are x\X2 and X2X4, 
respectively. The permutation Kf is in both cases the identity and with |/| = 6, |7| = 7 we have 
an overall minus sign 

= _i. (47) 

Therefore, the right hand side of eq. ((421 equals — x\X2 —X2X4, showing the agreement with the 
result of eq. (|46|) . 

Eq. (1371) together with eq. (|34l) allows to determine the first Symanzik polynomial U from 
the Laplacian of the graph. We may ask if also the polynomial f$ can be obtained in a similar 
way. We consider again a graph G with n internal edges [e\ , e n ), r internal vertices (vi , . . . , v r ) 
and m external legs. We attach m additional vertices v r +i, v r + m to the ends of the external legs 
and view the former external legs as additional edges e n+ \, e n+m . This defines a new graph 
G. As before we associate the parameters x\ to the edges e l (\ <i< n) and new parameters Zj to 
the edges e n+ j (1 < j < m). The Laplacian of G is a (r + m) x (r+m) matrix. Now we consider 
the polynomial 

W(xi,...,x n ,zi,...,z m ) = detL(G) [r+l,...,r + m]. (48) 

W is a polynomial of degree r = n — I + I in the variables x\ and Zj. We can expand W in 
polynomials homogeneous in the variables zf 

W = y l /(°) + w W + W {i) + ... + W {m \ (49) 
where W W is homogeneous of degree k in the variables Zj- We further write 

w (k) = £ ^l^{ Xll ... 1 x n ) ZJl .^ Jk . (50) 

Uu-Jk) 
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The sum is over all indices with 1 < j\ < ... < jj < m. The j \ are homogeneous polyno- 

mials of degree r — kin the variables x/. For W (°) and one finds 

m 

w^ = o, w^ = 9c(x h ...,x n ) J^zj, (5i) 

7=1 

therefore 

W = xi...x„^ } ^,...,^, (52) 
for any j e {1, ...,m}. Jo is related to TV ^ 2 ^: 




The proof of eqs. (|5TI)- (|5"3~1) follows from the all-minors matrix-tree theorem. The all-minors 
matrix-tree theorem states 

W(xi,...,x n ,zi,...,z m ) = W a p ( 54 ) 

with / = (r+ 1 , . . . , r+m) and aj = xj if ej is an internal edge or aj = Zj- n if ej is an external edge. 
The sum is over all m-forests of G, such that each tree in an m-forest contains exactly one of the 
added vertices v r +i, v r+m . Each m-forest has m connected components. The polynomial W (°) 
by definition does not contain any variable Zj- TV ^ would therefore correspond to forests where 
all edges connecting the external vertices e r+ i, e r+m have been cut. The external vertices 
appear therefore as isolated vertices in the forest. Such a forest must necessarily have more 
than m connected components. This is a contradiction with the requirement of having exactly m 
connected components and therefore < wM°) = 0. Next, we consider Each term is linear 

in the variables zj. Therefore (m — 1) vertices of the added vertices v r+ i, v r+m appear as 
isolated vertices in the m-forest. The remaining added vertex is connected to a spanning tree of 
G. Summing over all possibilities one sees that W W is given by the product of (zi + ... +z m ) 
with the Kirchhoff polynomial of G. Finally we consider W ^ . Here, (m — 2) of the added 
vertices appear as isolated vertices. The remaining two are connected to a spanning 2-forest of 
the graph G, one to each tree of the 2-forest. Summing over all possibilities one obtains eq. (|53l . 

Eq. (1521 and eq. (1531 together with eq. (|28l ) allow the computation of the first and second 
Symanzik polynomial from the Laplacian of the graph G. This provides a third method for the 
computation of the graph polynomials U and f . 

As an example we consider the double-box graph of fig. [T] We attach an additional vertex to 
every external line. Fig. |9] shows the labelling of the vertices and the Feynman parameters for the 
graph G. The Laplacian of G is a 10 x 10-matrix. We are interested in the minor, where - with 
the labelling of fig. |9]- we delete the rows and columns 7, 8, 9 and 10. The determinant of this 



16 



V 8. 



3 



X2 


x 5 




















X-j 




x 3 




v 3 


^6 l>5. 



v 7 



'Z3 



Figure 9: The labelling of the vertices and the Feynman parameters for the "double box"-graph. 

minor reads 

W =det L[7,8,9,10] = 



X\ +X 3 +Zi 


-xi 


~X3 











-X\ 


X\ +X 2 +Z2 





~X2 








-X3 





X3+X4+X 6 


—X4 


-x 6 










—X4 


X2+X4+X5 





-*5 








-x 6 





X(,+Xj +Z4 


-xi 













-xi 


X 5 +X7+Z3 




+ <W (3) + W 


(4)_ 









(55) 

For the polynomials W W and W ^ one finds 

= (21+22 + 23+24) 

(^1X2^3X5X6 + X{X2X3X5Xj + XiX2X3X(,Xj + X1X2X4X5X6 + 
J rX\X2X4X(,Xq + X^X^XfrX-j + X1X3X4X5X6 + X1JC3X4X5X7 + X\X3X4X(yX-] 
+X1X3X5X6X7 + X2X3X4X5X6 + X2X3X4X$X-] + X2X3X4X(,X-] + X2X3X5X6X7) , 
W (2) = {z\ J rZ4){z2 J cZ3)X2X3X5X< 3 

+ (Z1 +Z2) (Z4 + Z3) {x\X2X3X-] +X1X2X4X7 +XiX2X(,Xj + X1X3X4X7 + X1X3X5X7 + X1X4X5X6 
+X 1X4X5X7 +X 1X4X5X7 + X 1X5X5X7 +^2^3^4X7) 

+21 (Z2 +Z3 +Z4)X2 (X3X5X7 +X4X5X6 +X4X5X7 + X4X 6 X7 + X5X6X7) 
+22(21+23 +Z4)X3 (x2XeX7 + X4X5X5 + X4X5X7 + X4^7 + ^5^7) 
+Z3(zi+Z2+^4)^6 (X1X2X3 + X1X2X4+X1X3X4+X1X3X5 +Z2X3X4) 

+24(21 +22+23)^5 (X1X2X3 -\- X\X2X4 -\- X\X2X(j -\- X 1X3X4 -\- X2X3X4) . (56) 

With the help of eq. (1521 ) and eq. (|53l ) and using the kinematical specifications of eq. (fTOb we 
recover U and of eq. (PT4i 

We would like to make a few remarks: The polynomial W is obtained from the determinant 
of the matrix L = L(G) [r + 1 , . . . , r +m] . This matrix was constructed by attaching one additional 
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G G G 

Figure 10: The left picture shows a one-loop graph with three external legs. The middle pictures 
shows the associated graph G, where three additional vertices have been attached to the external 
legs. The right picture shows the graph G associated to G, where all external legs have been 
joined in one additional vertex. 

vertex to each external line. Then one deletes the rows and columns corresponding to the newly 
added vertices from the Laplacian of the graph obtained in this way. There are two alternative 
ways to arrive at the same matrix L: 

The first alternative consists in attaching to the original graph G a single new vertex v r +i, 
which connects to all external lines. This defines a new graph G, which by construction no 
longer has any external lines. As before we associate variables zi, z m to the edges connected 
to v r +i. Fig. [10] shows an example for the graphs G and G associated to a one-loop graph with 
three external legs. The Laplacian of G is a (r + 1) x (r + l)-matrix. It is easy to see that 

L = L(G)[r+l}. (57) 

From eq. (1371) we see that the polynomial W is nothing else than the Kirchhoff polynomial of 
the graph G: 

W(G) = 3c(G)=detL(G)[j], (58) 

where j is any number between 1 and r + 1 . 

For the second alternative one starts from the Laplacian of the original graph G. Let n be a 
permutation of (1, r). We consider the diagonal matrix diag (z n m, ...,z,m) . We can choose 
the permutation n such that 

= j, if the external line j is attached to vertex v,-. (59) 

We then have 

L = L(G) + dmg(z n ^),...,z n (r))\ Zm+l= _ =Zr=0 - (60) 

6 Deletion and contraction properties 

In this section we study two operations on a graph: the deletion of an edge and the contraction of 
an edge. This leads to a recursive algorithm for the calculation of the graph polynomials 11 and 
J . In addition we discuss the multivariate Tutte polynomial and Dodgson's identity. 
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Figure 1 1 : A one-particle-reducible graph: The edge e^ is called a bridge. Deleting 63 results in 
two connected components. 




Figure 12: The left picture shows the graph G of the two-loop two-point function. The middle 
picture shows the graph G — es, where edge es has been deleted. The right picture shows the 
graph G/es, where the two vertices connected to e^ have been joined and the edge e$ has been 
removed. 

In graph theory an edge is called a bridge, if the deletion of the edge increases the number of 
connected components. In the physics literature the term "one-particle-reducible" is used for a 
connected graph containing at least one bridge. The contrary is called "one-particle-irreducible", 
i.e. a connected graph containing no bridges. Fig.QT]shows an example. The edge is a bridge, 
while the edges e\, e^, e^ and e$ are not bridges. Note that all edges of a tree graph are bridges. 
An edge which is neither a bridge nor a self-loop is called a regular edge. For a graph G and a 
regular edge e we define 

G/e to be the graph obtained from G by contracting the regular edge e, 
G — e to be the graph obtained from G by deleting the regular edge e. (61) 

Fig. d2l shows an example. If the graph G has loop number / it follows that G — e has loop 
number (I — I), while G/e has loop number I. This follows easily from the formula I = n — r + 1 
for a connected graph: G — e has one edge less, but the same number of vertices, while G/e has 
one edge and one vertex less. 

Let us now study the behaviour of the Laplacian under these operations. Under deletion the 
Laplacian behaves as 

L(G-e k ) = L(G)\ Xk=0 , (62) 

i.e. the Laplacian of the graph G — e^ is obtained from the Laplacian of the graph G by setting the 
variable to zero. The behaviour of the Laplacian under contraction is slightly more compli- 
cated: As before we consider a graph with r vertices. Assume that edge e^ connects the vertices 
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v a and v r . The Laplacian L(G/e k ) is then a (r — 1) x (r — l)-matrix with entries 



L(G/e k ) 



I 



L(G) fla +L(G) rr + L(G) ar + L(G), a , if i = j = a, 

L{G)aj + L(G) rj , ifi = a,j^a, 

L(G) ia +L(G) ir , if j = a,i^ a, 

L{G)ij, otherwise. 



(63) 



Therefore the Laplacian of L(G/e k ) is identical to the minor L(G)[r] except for the row and 
column a. The choice that the edge e k is connected to the last vertex v r was merely made to 
keep the notation simple. If the edge connects the vertices v a and vt> with a < b one deletes from 
L(G) row and column b and modifies row and column a analogously to the formula above with 
b substituted for r. In particular we have B4l 



The deletion/contraction operations can be used for a recursive definition of the graph polynomi- 
als. For any regular edge e k we have 



The recursion terminates when all edges are either bridges or self-loops. This is then a graph, 
which can be obtained from a tree graph by attaching self-loops to some vertices. These graphs 
are called terminal forms. If a terminal form has r vertices and I (self-) loops, then there are 
(r — 1) "tree-like" propagators, where the momenta flowing through these propagators are linear 
combinations of the external momenta pi alone and independent of the independent loop mo- 
menta kj. The momenta of the remaining / propagators are on the other hand independent of the 
external momenta and can be taken as the independent loop momenta kj, j = 1 , . . . , /. Let us agree 
that we label the (r— 1) "tree-like" edges from 1 to r— 1, and the remaining / edges by r, n 
with n = r+ 1 — 1. We further denote the momentum squared flowing through edge j by qj. For 
a terminal form we have 



L(G/e k ) [a] 



L(G) [a,b]. 



(64) 



U{G) 



u(G/e k )+x k u(G-e k ), 
f (G/e k )+x k !Fo(G-e k ). 



(65) 




(66) 



In the special case that the terminal form is a tree graph, this reduces to 




(67) 



The Kirchhoff polynomial has for any regular edge the recursion relation 



x k 9C(G/e k ) + 9C(G-e k ). 



(68) 
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Note that the factory appears here in combination with the contracted graph G/e k . The recursion 
ends again on terminal forms. For these graphs we have with the conventions as above 

X =xi...x r -i. (69) 

The recursion relations eq. (1651) and eq. (1681) are proven with the help of the formulae, which 
express the polynomials 11 and %i in terms of spanning trees. For one uses the corresponding 
formula, which expresses this polynomial in terms of spanning 2-forests. As an example consider 
the polynomial 11 and the set of all spanning trees. This set can be divided into two sub-sets: 
the first sub-set is given by the spanning trees, which contain the edge e k , while the second 
subset is given by those which do not. The spanning trees in the first sub-set are in one-to- 
one correspondence with the spanning trees of G/e k , the relation is given by contracting and 
decontracting the edge e k . The second subset is identical to the set of spanning trees of G — e k . 
The graphs G and G — e k differ by the edge e k which has been deleted, therefore the explicit 
factor Xk in front of u(G — e k )- 

Eq. (|65l) and eq. (l66l) together with eq. ((281) provide a fourth method for the computation of 
the graph polynomials 11 and J . 

We now look at a generalisation of the Kirchhoff polynomial satisfying a recursion relation 
similar to eq. (|68l) . For a graph G - not necessarily connected - we denote by S the set of all 
spanning sub-graphs of G, i.e. sub-graphs H of G, which contain all vertices of G. It is not 
required that a spanning sub-graph is a forest or a tree. As before we associate to each edge e, 
a variable Xj. We will need one further formal variable q. We recall that the loop number of a 
graph G with n internal edges and r vertices is given by 

I = n-r + k, (70) 

where k is the number of connected components of G. We can extend the definition of the 
deletion and contraction properties to edges which are not regular. If the edge e is a self loop, 
we define the contracted graph G/e to be identical to G — e. The multivariate Tutte polynomial 
is defined by [45] 

Z(q,x h ...,x n ) = £<Z* (H) n^ ( 71 ) 

HeS eieH 

It is a polynomial in q and x\, x n . The multivariate Tutte polynomial generalises the standard 
Tutte polynomial I146I - I5T1I . which is a polynomial in two variables. For the multivariate Tutte 
polynomial we have the recursion relation 

Z(G) = x k z(G/e k ) + z(G-e k ), (72) 

where e k is any edge, not necessarily regular. The terminal forms are graphs which consists 
solely of vertices without any edges. For a graph with r vertices and no edges one has 

Z = q r . (73) 

The multivariate Tutte polynomial starts as a polynomial in q with q k if G is a graph with k 
connected components. If the graph G is not connected we write G = (Gi, G k ), where G\ to 
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Gk are the connected components. For a disconnected graph the multivariate Tutte polynomial 
factorises: 

z(G) = z{G l )...z(G k ). (74) 

Some examples for the multivariate Tutte polynomial are 

qx + q 2 , 
q(x+l), 

q (xix 2 + xi + x 2 ) + q 2 ■ (75) 

If G is a connected graph we recover the Kirchhoff polynomial %i (G) from the Tutte polynomial 
z(G) by first taking the coefficient of the linear term in q and then retaining only those terms 
with the lowest degree of homogeneity in the variables x, . Expressed in a formula we have 

9C(xi,...,x n ) = lim lim X l ~ r q~ l Z (q,Xxi,...,Xx n ) . (76) 

To prove this formula one first notices that the definition in eq. (fTTI) of the multivariate Tutte 
polynomial is equivalent to 

Z(q, Xl ,...,x n ) = q^q m ]\-. (77) 

One then obtains 

l l - r q- l z(qMu...Mn) = l^q^-^-^Uxi. (78) 

HeS e t eH 

The limits q — > and X — > select k[H) = 1 and 1(H) = 0, hence the sum over the spanning 
sub-graphs reduces to a sum over spanning trees and one recovers the Kirchhoff polynomial. 

At the end of this section we want to discuss Dodgson's identity^. Dodgson's identity states 
that for any n x n matrix A and integers i, j with 1 <i,j<n and i ^ j one has H52H53I 

det(A)det(A[/,j]) = det(A[/])det(A[j]) -det(A[/; j])dct(A[j;i]) ■ (79) 

We remind the reader that A [i, j] denotes a (n — 2) x (n — 2) matrix obtained from A by deleting 
the rows and columns i and j. On the other hand A[i; j] denotes a (n — 1) x (n — 1) matrix which 
is obtained from A by deleting the z'-th row and the j-th column. The identity in eq. (1791 ) has 
an interesting application towards graph polynomials: Let e a and be two regular edges of a 
graph G, which share one common vertex. Assume that the edge e a connects the vertices v,- 
and Vk, while the edge eb connects the vertices vj and v^. The condition i ^ j ensures that after 
contraction of one edge the other edge is still regular. (If we would allow i = j we have a multiple 

'Dodgson's even more famous literary work contains the novel "Alice in wonderland" which he wrote using the 
pseudonym Lewis Carroll. 



z(~) = 

Z(Q) = 

*(0) = 
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edge and the contraction of one edge leads to a self-loop for the other edge.) For the Kirchhoff 
polynomial of the graph G — e a — e b we have 

X(G-e a -e b ) = &QtL{G-e a -e b )[k\. (80) 

Let us now consider the Kirchhoff polynomials of the graphs G/e a — e b and G/eb — e a . One finds 

9c(G/e a -e b ) = detL(G-e a -e b )[i,k], 

<K(G/e h -e a ) = detL(G-e a -e b )[j,k). (81) 

Here we made use of the fact that the operations of contraction and deletion commute (i.e. 
G/e a — eb = (G — eb)/e a ) as well as of the fact that the variable x a occurs in the Laplacian 
of G only in rows and columns i and k, therefore L(G — e b ) [i, k] = L (G — e a — eb) [i, k) . Finally 
we consider the Kirchhoff polynomial of the graph G/e a /eb, for which one finds 

9C(G/e a /e b ) = detL(G-e a -e b )[i,j,k). (82) 

The Laplacian of any graph is a symmetric matrix, therefore 

detL(G-e a -e b )[i,k;j,k] = detL(G-e a -e b )[j,k;i,k]. (83) 

We can now apply Dodgson's identity to the matrix L(G — e a — e b )[k]. Using the fact that 
L (G — e a — e b ) [i, k; j, k] = L(G) [z, k; j, k] one finds Il54l 

% (G/e a - e b ) % (G/e h — e a ) — % (G — e a — e h ) K (G/e a /e b ) = (detL(G) [f , k; j, k] f . (84) 

This equation shows that the expression on the left-hand side factorises into a square. The ex- 
pression on the right-hand side can be re-expressed using the all-minors matrix-tree theorem as 
a sum over 2-forests, such that the vertex v& is contained in one tree of the 2-forest, while the 
vertices v\ and Vj are both contained in the other tree. 

Expressed in terms of the first Symanzik polynomial we have 

u{G/e a -e b )u{G/e b -e a )-u{G-e a -e b )u{G/e a /e b ) = ( —) ■ (85) 
The expression Ai is given by 

Feq .(i,kUm e t £F 

The sum is over all 2-forests F = (T\,T2) of G such that Vi,Vj G T\ and v\ E T%. Note that each 
term of Ai contains x a and x b . The factorisation property of eq. (|85l) plays a crucial role in the 
algorithms of refs. (6]-[8]|. 

A factorisation formula similar to eq. (|85l ) can be derived for an expression containing both 
the first Symanzik polynomial U and the polynomial Jq. As before we assume that e a and e b 
are two regular edges of a graph G, which share one common vertex. The derivation uses the 
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results of sect. [5] and starts from eq. (f8~4l) for the graph G associated to G. Eq. (1581) relates then 
the Kirchhoff polynomial of G to the W -polynomial of G. The W -polynomial is then expanded 
in powers of z. The lowest order terms reproduce eq. (|85l) . The next order yields 

11 (G/e a - e b ) ?q (G/e b - e a ) -u(G-e a - e b ) Jo (G/e a /e b ) (87) 

A 2 " 



+ J (G/e a - e b ) 11 (G/e b - e a ) - Jo (G - e a - e b ) 11 (G/e a /e b ) 

The quantity A 2 appearing on the right-hand side is obtained from the all-minors matrix-tree 
theorem. We can express this quantity in terms of spanning three-forests of G as follows: Let 
us denote by q , ^ l, ^ , ' ,k ^ the set of spanning three-forests (Ti,Ti, T$) of G such that v/,V/ G T\ and 

Vk G ?3. Similar we denote by T 3 ^ J '^ the set of spanning three-forests (T\,T2, T3) of G such that 
V; G 7\ , Vj G T2 and G ?3 . Then 

* = e e (^) n * 

(T h T2,T3)e<T 3 {i > J > k) v c£ T hV d €T 2 \ " / e t $(T u T 2 ,T 3 ) 

e e f^#) n *• (88) 

(r 1 ,r 2) r 3 )e<r 3 ((y) '' , '* ) v ^er 2 V j" / e ,i{T h T 2 ,T 3 ) 

In this formula we used the convention that the momentum /?j equals zero if no external leg is 
attached to vertex vj. Expanding the W -polynomial to order z 4 we have terms of order z 2 squared 
as well as terms which are products of order z with order z 3 . We are interested in an expression 
which arises from terms of order z 2 squared alone. In this case we obtain a factorisation formula 
only for special kinematical configurations. If for all external momenta one has 

(Ph ■ Ph) ■ {Ph ■ Pi 4 ) = {Ph-Ph)-{Ph-Pi4)i hih, hih G { 1, m} (89) 

then 

fo(G/e a -e b )Fo(G/e b -ea)-!F (G-e a -e b )!F {G/e a /e b ) = ( . (90) 



Eq. ([891 is satisfied for example if all external momenta are collinear. A second example is given 
by a three-point function. In the kinematical configuration where 

{p 2 l ) 2 +{p 2 2 ) 2 +{p 2 3 ) 2 -2p 2 p 2 -2p 2 p 2 -2p 2 p 2 = 0, (91) 

eq. (l89l) is satisfied. 



7 Duality 

We have seen that the Kirchhoff polynomial and the first Symanzik polynomial 11 of a graph 
G with n edges are related by the equations (|34~l) : 

ll(xi,...,x n ) =x\...x n 3C (—,...,— ) , %:(xi,...,x n ) =Xi...X„ 11 (—,...,— 

\X\ X n J \X\ %n 



24 




Figure 13: The "crossed double-box"-graph can be drawn as a planar graph. 




^5 ^3,3 

Figure 14: The 'smallest' non-planar graphs. 

In this section we will ask if one can find a graph G* with n edges such that %_ (G*) = U (G) and 
(G) = U (G*). Such a graph G* will be called a dual graph of G. In this section we will show 
that for a planar graph one can always construct a dual graph. The dual graph of G need not be 
unique, there might be several topologically distinct graphs G* fulfilling the above mentioned 
relation. In other words two topologically distinct graphs G\ and G2 both of them with n edges 
can have the same Kirchhoff polynomial. 

In this section we consider only graphs with no external lines. If a graph contains an external 
line, we can convert it to a graph with no external lines by attaching an additional vertex to the 
end of the external line. This is the same operation as the one considered already in fig. [9j A 
graph is called planar if it can be embedded in a plane without crossings of edges. We would like 
to note that the "crossed double-box"-graph shown in fig. [T3]is a planar graph. The right picture 
of fig. [T3l shows how this graph can be drawn in the plane without any crossing of edges. 

Fig .[141 shows two examples of non-planar graphs. The first graph is the complete graph with 
five vertices £5. The second example is denoted £3,3. A theorem states that a graph G is planar 
if and only if none of the graphs obtained from G by a (possibly empty) sequence of contractions 
of edges contains ^5 or ^3,3 as a sub-graph ll55l - [57l . 

Each planar graph G has a dual graph G* which can be obtained as follows: 

• Draw the graph G in a plane, such that no edges intersect. In this way, the graph divides 
the plane into open subsets, called faces. 

• Draw a vertex inside each face. These are the vertices of G*. 
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Figure 16: An example showing that different embedding s of a planar graph G into the plane 
yield different dual graphs G\ and G\. 

• For each edge e\ of G draw a new edge e* connecting the two vertices of the faces, which 
are separated by e,-. The new edges e* are the edges of G*. 

An example for this construction is shown in fig. \\5\ We note from the construction of the dual 
graph, that for each external edge in G there is a self-loop in G* and that for each self-loop in G 
there is an external edge in G*. 

If we now associate the variable Xi to the edge e\ of G as well as to the edge e* of G* we have 

K(G*) = U(G), K(G) = U(G*). (92) 

It is important to note that the above construction of the dual graph G* depends on the way, how 
G is drawn in the plane. A given graph G can have several topologically distinct dual graphs. 
These dual graphs have the same Kirchhoff polynomial. An example is shown in fig. [16] For this 
example one finds 

Ki{G) = U{G\) = U{G\) = (x\+X2)(x3+X4){X5+X 6 )+X3X4(x5+Xe) + (x3+X4)X5X6, 

11(G) = 3C(G|) = K(G* 2 )=x l X2(x 3 +X4+x 5 +x 6 ) + (x l +X2)(x3+x 4 )(x 5 +x 6 ). (93) 
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Figure 17: A graph G. 



8 Matroids 



In this section we introduce the basic terminology of matroid theory. We are in particular in- 
terested in cycle matroids. A cycle matroid can be constructed from a graph and it contains the 
information which is needed for the construction of the Kirchhoff polynomial and therefore as 
well for the construction of the first Symanzik polynomial 11 . In terms of matroids we want to 
discuss the fact, that two different graphs can have the same Kirchhoff polynomial. We have 
already encountered an example in fig. [16] We review a theorem on matroids which determines 
the classes of graphs whose Kirchhoff polynomials are the same. For a detailed introduction to 
matroid theory we refer to H581I59I . 

We introduce cycle matroids by an example and consider the graph G of fig. [FT] The graph G 
has three vertices V = {vi, V2, V3} and four edges E = {ei, e% ^3, e<\}. The graph has five span- 
ning trees given by the sets of edges {ei, e^}, {e\, e<x}, {^2, £3}, {^2, £4], {^3, £4], respectively. 
We obtain the Kirchhoff polynomial = x\X'$-\-x\X4 +^2^3 +-*2-* : 4 +-£3*4. The incidence matrix 
of a graph G with r vertices and n edges is a r x ^-matrix B — (by), defined by 

J 1, if ej is incident to v,-, , Q .. 
bij ~ \ 0, else. (y4j 

For the graph G of fig. [FT] the incidence matrix reads 





e 2 


<?3 


<?4 


1 


1 


1 










1 


:) 


1 


1 








where we indicated that each column vector corresponds to one edge of the graph. Let us focus 
on the set of these four column vectors. We want to consider all subsets of these vectors which are 
linearly independent over Z2. Obviously the set of all four vectors and the set of any three of the 
given vectors are linearly dependent over Z2. Furthermore the first two columns corresponding 
to {e\,e{\ are equal and therefore linearly dependent. Hence the linearly independent subsets 
are all sets with only one vector and all sets with two vectors, except for the just mentioned one 
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consisting of the first and second column. For each set of linearly independent vectors let us now 
write the set of the corresponding edges. The set of all these sets shall be denoted X. We obtain 



{0, {ei} , {e 2 } , {e 3 } , M , 

{<?!, e 3 }, j>l, <?4}, j>2, <? 3 }, j>2, ^4}, j>3, e 4 }}. 



(95) 



The empty set is said to be independent and is included here by definition. Let us make the 
important observation, that the sets in X which have two elements, i.e. the maximal number of 
elements, are exactly the sets of edges of the spanning trees of the graph given above. 

The pair {E, X) consisting of the set of edges E and the set of linearly independent sets X is 
an example of a matroid. Matroids are defined as ordered pairs (E, X) where E is a finite set, the 
ground set, and where Xis a collection of subsets of E, called the independent sets, fulfilling the 
following conditions: 

1. 0GX. 

2. If / G X and /' C /, then /' G X 

3 . If h and I2 are in X and |/i | < \h\, then there is an element eofh— h such that I\ U e G X. 

All subsets of E which do not belong to X are called dependent. The definition goes back to 
Whitney who wanted to describe the properties of linearly independent sets in an abstract way. 
In a similar way as a topology on a space is given by the distinction between open and closed sets, 
a matroid is given by deciding, which of the subsets of a ground set E shall be called independent 
and which dependent. A matroid can be defined on any kind of ground set, but if we choose E 
to be a set of vectors, we can see that the conditions for the independent sets match with the 
well-known linear independence of vectors. 

Let us go through the three conditions. The first condition simply says that the empty set 
shall be called independent. The second condition states that a subset of an independent set is 
again an independent set. This is fulfilled for sets of linearly independent vectors as we already 
have seen in the above example. The third condition is called the independence augmentation 
axiom and may be clarified by an example. Let the sets 



h 










1 





' 


w 


w 



h = 



/1\ 






1 
1 

V 1 / 






V 1 / 



be linearly independent sets of vectors, where I2 has more elements than I\ . I% — l\ is the set of 
vectors in I2 which do not belong to I\. The set I2 — h contains for example e = (1, 0, 0, l) r 
and if we include this vector in I\ then we obtain again a linearly independent set. The third 
condition states that such an e can be found for any two independent sets with different numbers 
of elements. 

The most important origins of examples of matroids are linear algebra and graph theory. The 
cycle matroid (or polygon matroid) of a graph G is the matroid whose ground set E is given by 
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the edges of G and whose independent sets X are given by the linearly independent subsets over 
1a2 of column vectors in the incidence matrix of G. We can convince ourselves, that X fulfills the 
conditions layed out above. 

Let us consider the bases or maximal independent sets of a matroid (E, X). These are the sets 
in X which are not proper subsets of any other sets in X The set of these bases of a matroid shall 
be denoted CB and it can be defined by the following conditions: 

1 . <B is non-empty. 

2. If B[ and B2 are members of (B and x E B\ — B2, then there is an element y of B2 — #i such 
that (Bi -x) Uy E <B. 

One can show, that all sets in <B have the same number of elements. Furthermore Xis uniquely 
determined by (B : it contains the empty set and all subsets of members of (2 . 

Let M = (E, X) be the cycle matroid of a connected graph G and let <B (M) be the set of bases. 
Then one can show that <B (M) consists of the sets of edges of the spanning trees in G. In other 
words, T is a spanning tree of G if and only if its set of edges are a basis in (B(M). We can 
therefore relate the Kirchhoff polynomial to the bases of the cycle matroid: 

sc = e n*< ( % ) 

Bj£-B{M) e/eBj 

The Kirchhoff polynomial of G is called a basis generating polynomial of the matroid M associ- 
ated to G. The Kirchhoff polynomial allows us to read off the set of bases $ (M) . Therefore two 
graphs without any self-loops have the same Kirchhoff polynomial if and only if they have the 
same cycle matroid associated to them. 

Let us cure a certain ambiguity which is still left when we consider cycle matroids and Kirch- 
hoff polynomials and which comes from the freedom in giving names to the edges of the graph. 
In the graph of fig. [17] we could obviously choose different names for the edges and their Feyn- 
man parameters, for example the edge e2 could instead be named and vice versa. As a conse- 
quence we would obtain a different cycle matroid where compared to the above one, takes the 
place of e2 and vice versa. Similarly, we would obtain a different Kirchhoff polynomial, where 
the variables X2 and X3 are exchanged. Of course we are not interested in any such different 
cases which simply result from a change of variable names. Therefore it makes sense to consider 
classes of isomorphic matroids and Kirchhoff polynomials. 

Let Mi and A/2 be two matroids and let E{M\) and EiMi) be their ground sets, respectively. 
The matroids M\ and M2 are called isomorphic if there is a bijection \\f from E{M\) to EiMi) 
such that is an independent set in M2 if and only if / is an independent set in M\ . The 
mentioned interchange of e2 and 63 in the above example would be such a bijection: V|/(e2) = 
\|/(c3) = e2, = £i for i = 1,4. The new matroid obtained this way is isomorphic to the 

above one and its independent sets are given by interchanging e2 and ej, in Xof eq. (|95l ). In the 
same sense we want to say that two Kirchhoff polynomials are isomorphic if they are equal up 
to bijections on their sets of variables. 

Now let us come to the question when the Kirchhoff polynomials of two different graphs 
are isomorphic, which means that after an appropriate change of variable names they are equal. 
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G' 



Figure 18: Vertex identification and cleaving 

From the above discussion and eq. (l96l it is now clear, that these are exactly the cases where the 
cycle matroids of the graphs are isomorphic. The question when two graphs have isomorphic 
cycle matroids was answered in the following theorem of Whitney If60l (also see ItoTlO which 
was one of the foundational results of matroid theory: 

Let G and H be graphs having no isolated vertices. Then the cycle matroids M(G) and 
M(H) are isomorphic if and only if G is obtained from H after a sequence of the following three 
transformations: 

1. Vertex identification: Let u and v be vertices of distinct components of a graph G. Then a 
new graph G is obtained from the identification of u and v as a new vertex w in G (see the 
transition from G to G in fig. [TBI . 

2. Vertex cleaving: Vertex cleaving is the reverse operation of vertex identification, such that 
from cleaving at vertex w in G we obtain u and v in distinct components of G (see the 
transition from G to G in fig. [T8l . 

3. Twisting: Let G be a graph which is obtained from two disjoint graphs G\ and G2 by 
identifying the vertices u\ of G\ and U2 of G2 as a vertex u of G and by identifying the 
vertices vi of G\ and V2 of G2 as a vertex v of G. Then the graph G is called the twisting 
of G about {w, v} if it is obtained from G\ and G% by identifying instead u\ with V2 and vi 
with U2 (see fig.[T9l). 

Proofs can be obtained from H58U60J - I621 . As a consequence of Whitney's theorem, the Kirchhoff 
polynomials %i (G) and (H) of the connected graphs G and H, both without any self-loops, are 
isomorphic if and only if G is obtained from H by a sequence of the above three transformations. 
For the transformations of vertex identification and vertex cleaving this is obvious from a well- 
known observation: If two distinct components G\ and G2 are obtained from G after vertex 
cleaving, then 9C (G) is the product of %i (G\) and (G2). Therefore any other graph G obtained 
from G\ and G2 after vertex identification has the same Kirchhoff polynomial (G ; ) = %i (Gi) ■ 
%. (G2) = K (G). The non-trivial part of the statement on the Kirchhoff polynomials of G and H 
is the relevance of the operation of twisting. In the initial example of fig. [16] the two graphs G\ 
and G* 2 can be obtained from each other by twisting. 
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Figure 19: Twisting about u and v 

There is an alternative definition of a matroid: Instead of specifying the set X of independent 
sets a matroid can be defined by a rank function, which associates a non-negative integer to every 
sub-set of the ground set. The rank function has to satisfy for all 5,5' C E the following three 
conditions: 

1. rk(S) < \S\. 

2. S' c S implies rk(S') < rk(S). 

3. rk(SUS')+ rk (5n5') <rk(S) + rk(S'). 

The independent sets are exactly those for which rk(5) = \S\ holds. For the cycle matroid of a 
graph G we can associate to a subset S of E the spanning sub-graph H of G obtained by taking 
all the vertices of G, but just the edges which are in S. In this case the rank of S equals the 
number of vertices of H minus the number of connected components of H. The multivariate 
Tutte polynomial for a matroid is defined by 

..,*„) = £<r rk(5) n^-- ( 97 ) 

SCE e,eS 

It is a polynomial in l/q and x\, x n . Since a matroid can be defined by giving the rank for 
each subset S of E, it is clear that the multivariate Tutte polynomial encodes all information of a 
matroid. For the cycle matroid of a graph G with r vertices the multivariate Tutte polynomial Z 
of the matroid is related to the multivariate Tutte polynomial Z of the graph by 

Z(q,xi,...,x n ) = q~ r z(q,xi,...,x n ). (98) 

For a matroid there are as well the notions of duality, deletion and contraction. Let us start with 
the definition of the dual of a matroid. We consider a matroid M with the ground set E and the 
set of bases $(M) = ,52, B n }. The dual matroid M* of M is the matroid with ground set 
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E and whose set of bases is given by <B (M*) = {E — B\, E — B2, ...,E — B n }. It should be noted 
that in contrast to graphs the dual matroid can be constructed for any arbitrary matroid. 

Deletion and contraction for matroids are defined as follows: Let us consider a matroid M = 
(E, X) with ground set E and Xbeing the set of independent sets. Let us divide the ground set E 
into two disjoint sets X and Y: 

E=XUY, lny = 0. (99) 

We denote by ly the elements of X, which are sub-sets of Y. The matroid M — X is then defined 
as the matroid with ground set E — X = Y and whose set of independent sets is given by Xy. 
We say that the matroid M — X is obtained from the matroid M by deleting X. The contracted 
matroid M/X is defined as follows: 

M/X = (M*-X)*, (100) 

i.e. one deletes from the dual M* the set X and takes then the dual. The contracted matroid M/X 
has the ground set E —X. With these definitions of deletion and contraction we can now state 
the recursion relation for the multivariate Tutte polynomial of a matroid: We have to distinguish 
two cases. If the one-element set {e} is of rank zero (corresponding to a self-loop in a graph) we 
have 

Z(M) = z(M-{e})+x e Z,(M/{e}). (101) 

Otherwise we have 

z(M) = z(M-{e}) + ^Z,(M/{e}). (102) 

The recursion terminates for a matroid with an empty ground set, in this case we have Z = 1 . 
The fact that one has a different recursion relation for the case where the one-element set {e} is 
of rank zero is easily understood from the definition of Z and the relation to graphs: For a cycle 
matroid Z differs from Z by the extra factor q~ r , where r is the number of vertices of the graph. 
If e is a self-loop of G, the contracted graph G/e equals G — e and in particular it has the same 
number of vertices as G. In all other cases the contracted graph G/e has one vertex less than G, 
which is reflected by the factor l/q in eq. (1 1 02b . 

9 Conclusions 

In this review we considered polynomials associated to Feynman graphs. The integrand of a 
Feynman loop integral is directly related to two graph polynomials, which are called the first and 
the second Symanzik polynomial. Upon a simple substitution the first Symanzik polynomial is 
related to the Kirchhoff polynomial of the graph. The graph polynomials have many interesting 
properties. We discussed the interpretation in terms of spanning trees and spanning forests, the 
all-minors matrix-tree theorem as well as recursion relations due to contraction and deletion of 
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edges. These properties in turn provide a variety of methods to compute the graph polynomials. 
We provided a detailed discussion of Dodgson's identity. This identity provides a factorisation 
formula for a certain combination of graph polynomials and can be used in practical calculations. 
In addition we discussed the generalisation of the two graph polynomials towards the multivariate 
Tutte polynomial and matroids. Matroid theory gives an answer to the question under which 
conditions two inequivalent graphs have the same Kirchhoff polynomial. 
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